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CN . Abstract 
bJOl We discuss our recent results on the statistical properties of damped 

^ ■ Lyman alpha systems (DLAs) at low redshift {z < 1.65) (Rao & Turnshek 

I 2000). Contrary to expectations, we found that the cosmological neutral 

QQ ■ gas mass density as traced by DLAs, ^dla, does not evolve from redshifts 

z ~ 4 to z ~ 0.5 and that extrapolation to z = results in a value that 
is a factor of ~ 6.5 times higher than what is derived from galaxies at 
^ , the current epoch using HI 21 cm emission measurements. We review 

the current status of HI measurements at low redshift and at the current 
epoch, and discuss possible causes of this discrepancy. 
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1. Introduction 

Q . It was recognized in the mid-eighties that damped Lyman alpha systems (DLAs) 

^ I reveal unique information about the evolutionary history of the Universe. Opti- 

^ . cal surveys of QSOs showed that DLAs trace the bulk of the observable neutral 

gas mass at high redshift, and therefore, that they can be used to study the for- 
mation and evolution of galaxies. A particularly compelling result was that the 
comoving neutral gas mass density in DLAs at redshift z ~ 3.5 was comparable 
to the luminous mass density observed in galaxies at the current epoch. There 
\ also seemed to be evidence for a decline in the neutral gas mass density from 

z « 3.5 to z ~ 1.7 that extrapolated to the value inferred from gas-rich galaxies 
at z = (Wolfe et al. 1986; Turnshek et al. 1989; Lanzetta et al. 1991; Rao 
&: Briggs 1993). This was interpreted as evidence for gas depletion due to star 
formation, and DLAs were thought to trace the evolution of the progenitors of 
present-day disk galaxies. 

However, several studies carried out during the second half of the last decade 
have shown that the interpretation is not so straightforward. Groundbased 
surveys can search for DLA lines only beyond z = 1.65 when the Lya line 
moves into the optical region of the electromagnetic spectrum. Thus, until the 
recent lUE and HST surveys for DLAs (Lanzetta, Wolfe, & Turnshek 1995; 
Rao, Turnshek, & Briggs 1995; Jannuzi et al. 1998; Rao k Turnshek 2000, 
henceforth, RT2000), about 70% of the Lya Universe was unexplored. The 
RT2000 survey, which we discuss in this contribution, more than doubled the 
number of known low-redshift DLAs. We used a bootstrapping method which 
relied on the statistics of Mgll absorption lines to identify DLAs. The results, 
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although still limited by small number statistics, show that ^dla is consistent 
with remaining constant from redshifts z 4 to z ft! 0.5 and that extrapolation 
to 2; = results in a value that is a factor of ~ 6.5 times higher than what is 
derived from galaxies at the current epoch. The DLA redshift number density is 
only marginally consistent with evolution between redshifts z ^ 4 and z « 0.5. 
This implies that DLAs may be a slowly evolving population that do not undergo 
significant gas depletion and, by extension, star formation. 

Apart from the low-redshift DLA statistics, results from metallicity studies, 
kinematics, and direct imaging of the DLAs have shown that the "DLA - disk 
galaxy" paradigm is inadequate. Pettini and collaborators (e.g. Pettini et al. 
1997, Pettini ct al. 1999, Pettini et al. 2000) have shown that there is no 
obvious tendency for DLA metallicities to increase with decreasing redshift and 
approach the solar value. However, there is considerable scatter in the individual 
metallicities which is most likely caused by the wide range of formation histories 
and galaxy types responsible for the DLAs. In fact, Pettini et al. (1999) have 
also concluded that the known DLAs are unlikely to trace the galaxy population 
responsible for the bulk of the star formation. 

Prochaska & Wolfe (1997; 1998) have used high-resolution spectra of the 
metal-line profiles associated with DLAs to argue that the kinematic profiles 
of the DLA absorbing regions are more consistent with models of rotating HI 
disks than any other single type of model. The kinematics have also been shown 
to be consistent with gas in-fall due to merging (Haehnelt, Steinmetz, & Ranch 

1998) , randomly moving clouds in a spherical halo (McDonald & Miralda-EscTidc 

1999) , and multiple gas disks in a common halo (Mailer et al. 2000). Therefore, 
it seems most likely that, consistent with other current findings, the kinematics 
of DLA absorbing regions arise from a mix of kinematic structures, including 
some rotating gaseous disks. 

There is now direct evidence from imaging studies that the morphological 
types of DLA galaxies are indeed mixed. Sec Rao &; Turnshek (1998), Turnshek 
et al. (2001), and Nestor et al. (2001, these proceedings) for ground-based 
imaging results on some of our low-redshift DLA fields. HST images of other low- 
redshift DLA absorbers have also revealed a mixed population (Le Brun et al. 
1997; Bouche et al. 2001). Turnshek, Rao, & Nestor (2001, these proceedings) 
compare the properties of low-redshift DLA galaxies to the properties of gas-rich 
galaxies at z = 0. 

Thus, there has been a shift in our interpretation of the nature of the DLA 
galaxy population, consistent with the idea that DLAs arise in giant hydrogen 
clouds that could be associated with any type of galaxy or protogalaxy (Kher- 
sonsky & Turnshek 1996). Here, we elaborate on some of the statistical results 
from our DLA survey. 



2. The Low-Redshift DLA Survey 

We recently completed an efficient non-traditional (but unbiased) HST survey to 
discover DLAs at redshifts z < 1.65 (RT2000). Our survey relied on observations 
of LycK absorption in identified Mgll systems (which can be identified from the 
ground for z > 0.1), where the incidence of Mgll is known as a function of the 
MgIIA2796 rest equivalent width, ^5^2796 g gteidel k Sargent 1992). The 
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Figure 1. A plot of Mgll Wq^'^^'^ versus Fell Wo^'^°°. The DLAs from the 
RT2000 sample are marked with open squares. The open circles represent previ- 
ously known 21 cm absorbers that were excluded from our unbiased Mgll sample. 
Left pointing arrows indicate upper limits to the measured value of Wq . The 
horizontal and vertical dashed lines identify the region for which Wo^''^^> 0.5 A 
and Wo^2«°° > 0.5 A; half of these are DLAs. 

empirical fact that all DLA absorbers have Mgll absorption (Turnshek et al. 
1989; Lu ct al. 1993; Wolfe et al. 1993; Lu & Wolfe 1994) greatly improves 
the efficiency of searches for DLA at low redshift. With our survey technique, 
we uncovered 12 DLA lines in 87 Mgll systems with Wq'^'^^^ > 0.3 A (a success 
rate of ~14%). Two more DLAs were discovered scrcndipitously. In total our 
survey increased the number of known low-redshift DLA absorbers more than 
two-fold (previously known low-redshift DLAs were mostly identified via 21 cm 
absorption). 

2.1. The Statistical Properties of Low- Redshift DLAs 

Since we observed the fraction of Mgll systems with DLA, we determined the 
incidence of DLAs at low redshift by bootstrapping from the Mgll statistics. By 
fitting Voigt damping profiles to the Lyct lines in our UV spectra, we deduced 
Nhi for each DLA system and determined the low-redshift cosmological neutral 
gas mass density of DLA absorbers. Specifically, our survey results indicate: 

(1) Approximately 50% of the systems with VFg^^Tge >o.5 A andFell Wg^^eoo 
>0.5 A have DLA absorption (Figure 1). Thus, we uncovered a new selection 
criterion for the identification of DLAs. Moreover, all the non-DLAs in this 
regime have HI column densities Nhi > 10^^ atoms cm~^. We also found that 
with the exception of the known z = 0.692 21 cm absorber towards 3C 286 
which has unusually low metal-line equivalent widths (T''l''f^^^^*' = 0.39 A and 
^A2600 ^ 0.22 A, Cohen et al. 1994), all of the DLAs lie in the upper-right 
region of the plot where M/g^^TQe > q.S A and Wo^2600 > q.S A. 

(2) The incidence of DLAs per unit redshift, no la, is observed to decrease 
with decreasing redshift (Figure 2). The observed trend in nuLA implies sig- 
nificant evolution only if the z = data point, which is derived from 21 cm 
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Figure 2. The number density redshift distribution of DLAs. Our low-redshift 
results are shown as solid squares (RT2000). A Malmquist bias correction factor 
has been applied to the high-redshift data of Wolfe et al. (1995) and the results are 
shown as solid circles. See RT2000 for details. The open circle is derived from the 
observed HI distribution in local spiral galaxies (Rao & Briggs 1993). The solid line 
hcis 7 = 2.5 and is forced to pass through the z = data point, the dashed line has 
7 = 1.5 and docs not include the z — data point. Extrapolation of this power-law 
to 2; = (the open triangle) results in a value that is ~ 3 times larger than the 
observed incidence at z = 0. 



observations of gas-rich spirals (Rao & Briggs 1993), is included in the analysis. 
Two power laws of the form noLAiz) = no(l + z)'^ are shown. For A = cos- 
mologies, riDLA is consistent with no intrinsic evolution if 7 = 1.0 for go = or 
if 7 = 0.5 for qq = 0.5. The solid line in Figure 2, which goes through the z = Q 
data point, has exponent 7 = 2.5 and indicates significant evolution. The dashed 
line, which excludes the z = data point, has 7 = 1.5 and is only marginally 
consistent with evolution. Extrapolation to 2 = (the open triangle) results in 
a value that is ~ 3 times larger than the observed incidence at 2; = 0. 

(3) The cosmological mass density of neutral gas in low-redshift DLA ab- 
sorbers, ^DLAi is found to be comparable to that observed at high redshift 
(Figure 3). The error bars are large because Q.j:,la{z) is very sensitive to the 
small number of systems which have the highest column densities. For go = 0-5, 
^dla{z) is a factor of ~ 6.5 times larger than the value for ilgas{z = 0) inferred 
from local 21 cm observations. 

(4) The HI column density distribution (CDD) of the low-redshift DLA ab- 
sorber population is very different in comparison to high-redshift DLA absorbers, 
and in comparison to the column density distribution inferred from local spirals 
(Figure 4). The low-redshift DLAs exhibit a significantly larger fraction of very 
high column density systems in comparison to determinations at both high red- 
shift and locally. At no redshift does the CDD fall-off in proportion to ^N^^. 
An ~iV^| fall-off is theoretically predicted for disk-like systems (e.g. Milgrom 
1988) and this is, in fact, observed locally in spiral samples (Rao & Briggs 1993; 
Zwaan et al. 1999). 
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Figure 3. The cosmological mass density of neutral gas as a function of redshift 
for Ho = 65 km Mpc"'^, go = 0.5, and A = 0. The data points for z > 1.5 
are from Storrie-Lombardi and Wolfe (2000). The open circle at « = is the local 
neutral gas mass density as measured by Rao & Briggs (1993) and the open square 
at X = is the local luminous mass density in stars (Fukugita, Hogan, & Peebles 
1998). 




Figure 4. The normalized cumulative column density distributions for the low- 
redshift DLA sample (RT2000), the high-redshift DLA sample (Wolfe et al. 1995), 
and local galaxies (Rao & Briggs 1993). 
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Despite the increase in the number of known low-redshift DLAs, the statis- 
tics are still dominated by small numbers leading to large error bars on noLA 
and i^DLA- It could be argued that these quantities are also consistent with a 
real decrease between z ~ 4 and z ~ 0.5. However, the CDD at low redshift is 
very significantly different from that at high redshift. There is only a 2.8% prob- 
ability that the two samples are drawn from the same parent population. Also, 
the low-redshift DLA sample is different from the local spiral galaxy population 
at the 99.99% confidence level. 



3. Discussion 

In summary, the observational results on DLAs over the past decade have shown 
that their redshift number density, cosmological mass density, and metallicities 
show little evolution from redshifts z 4to z 0.5. We have also seen that DLA 
galaxies cover a wide range of gas-rich morphological types from predominantly 
dwarf and low surface brightness galaxies to some spirals. 

That they are a slowly evolving population implies that either: (1) they 
do not trace the bulk of the galaxy population responsible for star formation, 
i.e., they are a different population than the galaxies that are used to determine 
the star formation history of the Universe, or (2) that the galaxy population 
responsible for star formation is a small subset of the DLA galaxy population. 
If the first possibility holds, this is probably the result of a selection effect in 
that optical surveys of QSOs preferentially probe lines of sight that are relatively 
dust free. These sight lines must avoid star forming regions that are enshrouded 
in metal-enriched dust (Pei, Fall, & Hauser 1999). Indeed, the models of Pei, 
Fall, & Hauser (1999) imply that as much as ~70% of the neutral gas mass is 
being missed by DLA surveys. However, the recent study of DLAs detected in 
a sample of radio-selected quasar spectra by Ellison et al. (2001) suggests that 
the effect might be less important. 

Our results on the incidence of DLAs show that udla is only marginally 
consistent with evolution from z ~ 4 to z ~ 0.5, and that the extrapolated value 
of nDLA{z) at z = is higher than that deduced from 21 cm emission data of 
nearby galaxies by a factor of ~ 3 (sec Figure 2). We also found that ^dla 
does not evolve from z ~ 4 to z ~ 0.5 and that the extrapolated value at z = 
is a factor of ~6.5 times higher than the neutral gas mass density deduced from 
gas-rich galaxies at the present epoch (Figure 3). Thus, the QSO absorption-line 
results arc highly inconsistent with the results from 21 cm observations of local 
galaxies. Recently, Churchill (2001) used HST archival spectra to determine the 
incidence of Mgll systems at z 0.05 and, based on the RT2000 method of 
bootstrapping from Mgll statistics and the assumption of no evolution in the 
Mgll-to-DLA statistics, derived a value for nDLA{z ~ 0.05) which is essentially 
equivalent to the RT2000 value for nDLA{z ~ 0.5). He further noted that if 
the CDD of DLAs did not evolve from z ~ 0.5 to z ~ 0.05, then Q.£)la{z ~ 
0.05) would also be consistent with the RT2000 result for ^dla{z « 0.5). The 
implication, then, would be that ^dla is constant from z 4 to nearly the 
current epoch, i.e., z ~ 0.05. If the Churchill (2001) results arc confirmed 
with follow-up observations of the corresponding Lyct lines then the discrepancy 
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between the QSO absorption line results and 21 cm emission measurements of 
local galaxies would have to be explained. 



3.1. At low redshift 

It is possible that we are so dominated by errors from the statistics of small 
numbers, that our sample just happened to have a higher fraction of the highest 
column density systems by sheer chance. As mentioned above, the Churchill 
(2001) result on ^dla{z = 0.05) was derived assuming that the CDD of DLAs 
did not evolve from z ~ 0.5 to z = 0.05. So that result is also subject to ques- 
tion in the same way. The only way around this is to conduct larger surveys 
for DLAs at low redshift. In addition to improving the statistical uncertainties, 
this would permit a better determination of the shape of the HI CDD at large 
column densities, effectively allowing a determination of the maximum HI col- 
umn density. Understanding this turn-down in column density is important for 
the determination of Q^LA- It is also worth reiterating that the redshift interval 
< z < 1.65 includes the most recent ~70% of the age of the Universe. We 
certainly require more than a dozen or so DLAs to do justice to this era. 

It might be argued that our technique of selecting from a Mgll sample 
biases the DLA sample towards higher column densities. However, we have 
shown that our data are not consistent with this (see RT2000). There is no trend 
between Mgll Wq'^'^^^ and HI column density. While we used a Wq'^'^^^ detection 
threshold of O.SA, all but one of our DLAs have Wq'^'^^^ > 0.6A (the one system 
has Wq*^''^^ = 0.52A). Thus, we believe that our empirically-determined cut-off 
in Wq'^'^^^ is sound. It was also suggested by Prank Briggs at this conference 
that a higher Wq^^^^ cut-off should be used to search for DLAs since they are 
probably more indicative of the presence of a DLA system. Moreover, the Steidel 
& Sargent (1992) Mgll study showed that systems with W(^2796 > ^ q ^ evolved 
away faster than those with Wq'^'^^^ > O.gA. This might, he argued, result in a 
smaller value for ^Idla{z « 0.5). We find that the results on ^dla{z ~ 0.5) 
from the two Wq^'^^^ threshold samples in our data set are consistent with each 
other; of course, the error bars for the Wq^^^^ > l.OA sample are even larger. 

It is also a concern that our QSO samples might be influenced by gravita- 
tional lensing bias; DLA galaxies magnify QSOs, preferentially introducing them 
into optically selected samples. Le Brun et al. (2000) have shown this effect to 
be less than 0.3 magnitudes for their sample of 7 DLA galaxies. They also point 
out that our QSO sample is brighter than theirs and that magnification bias 
may, therefore, be more important for our sample. The QSOs with the highest 
column density DLAs in our sample have not been imaged with HST, and so the 
possibility of lensing by the DLA galaxies has not been studied with the greatest 
possible sensitivity. On the other hand, our ground-based studies suggest that 
DLA galaxies are not highly luminous, and therefore, massive galaxies. 

In any case, it is clear that in addition to increasing DLA sample sizes, 
biases due to Mgll selection, dust obscuration, and the effect of gravitational 
lensing could be assessed more carefully. 
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3.2. At the present epoch 

We should also re-examine the z = statistics. In Rao Sz Briggs (1993) we used 
the best available optical luminosity functions known at that time in conjunction 
with empirical HI mass - optical luminosity relations of gas-rich galaxies to 
calculate the HI mass density locally. This is the value plotted in Figure 3. The 
upper error bar has been modified slightly as explained in Rao, Turnshek, &: 
Briggs (1995). Fall &: Pei (1993) derived approximately the same value using 
mean values of the local luminosity density and Mhi/Lb- This exercise was 
repeated by Natarajan & Pcttini (1997) using more recent optical luminosity 
functions and they confirmed the Rao k, Briggs (1993) result. For comparison, 
the results on ^gasiz = 0) were 2.4x lO-^/i"^ (Rao & Briggs 1993), 2.6x lO-^/i"^ 
(Fall & Pei 1993), and 2.5 x IQ-'^/i-i (Natarajan & Pettini 1997). 

Direct measurements of the local gas mass density using HI 21 cm emis- 
sion surveys have also been found to be consistent with the above results. The 
Arecibo HI Strip Survey (AHISS, Zwaan et al. 1997) performed the most sensi- 
tive search for local 21 cm emission to date, being 5 times more sensitive than the 
Arecibo Dual Beam Survey (ADBS) of Rosenberg &; Schneider (2000), although 
the latter covered a larger volume. The limiting column density of the Zwaan et 
al. survey was 10"*^^ cm^^ at the 5(T level for gas filling the telescope beam, and 
they had the capability of detecting HI masses of 6 x lO^/i~^M0 at 7h~^ Mpc. 
For comparison, HIPASS, which is a blind 21 cm survey of the southern sky, 
reached a limiting column density of 7 x 10^^ cm~^ at the 3(J level and a mass 
limit of ~ 7 X Wh'^ Mq at Ih'^ Mpc (Kilborn, Webster, k Staveley-Smith 
1999; Kilborn 2001). All three 21 cm surveys have been used to construct the 
local HI mass function. The faint end slope of the HI mass function has been 
the focus of much debate since it bears directly on the values of Q-gas{z = 0) and 
n{z = 0) (Figures 2 and 3). Moreover, a large faint end slope might imply the 
existence of a new dwarf galaxy population that has not been identified opti- 
cally. For comparison, Zwaan et al. (1997) derive a slope of a = 1.2, Rosenberg 
& Schneider (2001, these proceedings) get a = 1.5, while Kilborn (2001) derives 
a = 1.52 (where a is the standard Schechter luminosity ftmction parameter). 
The integral of the HI mass function is an estimate of the total HI mass density 
at ^; = 0, and its ratio with the critical mass density at the current epoch gives 
^Hi{z = 0). A correction for a neutral gas composition of 75% H and 25% He 
by mass then gives Qgasiz = 0). The results on 0,gas{z = 0) for the three sur- 
veys are 2.7 x IQ-'^h-^ (Zwaan et al. 1997), 3.8 x lO-'^/i^^ (Kilborn 2001), and 
4.4 X 10~^/i~^ (from data in Rosenberg k Schneider 2001, these proceedings). 

Thus, the value of 0,gas{z = 0) derived from optical luminosity functions is 
consistent with the value derived from HI mass functions. As expected, the HI 
mass functions with a ~ 1.5 result in slightly larger values for ^gas{z = 0), but 
only at the 2a level. The similarity between the HI mass function estimated 
from AHISS with that derived from optical luminosity functions (Briggs k Rao 
1993) led Zwaan et al. (1997) to conclude that there is no large population of 
Hl-rich galaxies (e.g. low surface brightness or dwarf) that have been missed by 
optical galaxy surveys. See also Zwaan, Briggs, k Sprayberry (2001) who show 
that the optical luminosity function of HI selected galaxies is in agreement with 
the luminosity function of optically selected late-type galaxies. However, based 
on the sensitivities of the three 21 cm surveys discussed above, the possibility of 
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there being objects with masses lower than ~ lO''' Mq that contain DLA column 
densities, but that have a mean column density less than r-^ 10^^ cm^^ within 
the telescope beam cannot be ruled out. 

Rosenberg &; Schneider (2001, these proceedings) have also suggested that 
the faint-end slope of the ADBS HI mass function, a = 1.52, is large enough 
to make the contribution of low-HI-mass galaxies important to the DLA cross- 
section at z = 0. They derive a value for n{z = 0) that is a factor of ~ 5.5 
higher than our estimate from a complete sample of 27 gas-rich galaxies (see Rao, 
Turnshek, & Briggs 1995 and Turnshek, Rao, & Nestor 2001, these proceedings). 
While the ADBS sample covers a larger range in HI masses (5 X IO^Mq to 
2 X IO^^Mq) compared to our sample (5 x lO'^M© to 2 x 10^°Mq), the cross- 
sectional areas of our galaxies at the 2 x 10^'^ atoms cm^^ HI column density 
contour are systematically lower by a factor of ~ 4. This might possibly account 
for the larger value of n{z = 0) that they derive from their data set. 



4. Concluding Remarks 

Observations over the past several years have led us to a new paradigm for 
DLAs, namely, that DLAs trace clouds of neutral gas whose properties have not 
changed significantly between z ~ 4 and z ~ 0.5, and that these clouds arise in a 
variety of gas-rich galaxy types. There exists a discrepancy of about a factor of 
6.5 between the cosmological neutral gas mass density at z = as inferred from 
DLA statistics and 21 cm emission measurements of local galaxies. Whether 
this discrepancy is due to small-number statistics of DLAs at low redshift, a real 
selection difference, or real evolution can only be resolved with more and better 
data at low redshift and at z = 0. 

This research has been supported by STScI grant GO-06577.01-95A and 
NASA LTSA grant NAG5-7930. 
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